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Abstract. An exact solution of Einstein's field equations in empty space first 
found in 1985 by Quevedo and Mashhoon is analyzed in detail. This solution 
generalizes Kerr spacetime to include the case of matter with arbitrary mass 
quadrupole moment and is specified by three parameters, the mass M, the angular 
momentum per unit mass a and the quadrupole parameter q. It reduces to the 
Kerr spacetime in the limiting case 5 = and to the Erez-Rosen spacetime when 
the specific angular momentum a vanishes. The geometrical properties of such 
a solution are investigated. Causality violations, directional singularities and 
repulsive effects occur in the region close to the source. Geodesic motion and 
accelerated motion are studied on the equatorial plane which, due to the reflection 
symmetry property of the solution, turns out to be also a geodesic plane. 



PACS number: 04.20.Cv 

1. Introduction 

The problem of describing the gravitational field of astrophysical bodies is of central 
importance in general relativity, both as an issue of principle and as a foundation for 
explaining the results of observations. It is an issue of principle because general 
relativity is believed to be the most accurate theory of the gravitational field. 
Consequently, Einstein's theory should accept the existence of exact solutions that 
correctly describe the gravitational field of realistic sources. On the other hand, the 
explanation of effects observed at the astrophysical level is extremely important. 

Astrophysical bodies are characterized in general by a non-spherically symmetric 
distribution of mass. In many cases, like ordinary planets and satellites, it is possible 
to neglect the deviations from spherical symmetry; it seems instead reasonable to 
expect that deviations should be taken into account in case of strong gravitational 
fields. 

The general metric describing the gravitational field of a rotating deformed mass 
was found in 1991 by Quevedo and Mashhoon [H [2] and involves an infinite set of 
gravitoelectric and gravitomagnetic multipoles. This is a stationary axisymmetric 
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solution of the vacuum Einstein's equations belonging to the class of Weyl-Lewis- 
Papapetrou [21 SI [S] and is characterized, in general, by the presence of a naked 
singularity. To capture the main properties of the general solution [U [5], we 
concentrate for the sake of simplicity on the special case of the general solution that 
involves only three parameters: the mass Af , the angular momentum per unit mass a 
and the mass quadrupole parameter q of the source. This special case was first found 
by Quevedo and Mashhoon in 1985 |61. Hereafter this solution will be denoted as the 
QM solution. 

The corresponding line element in prolate spheroidal coordinates {t,x,y^<j)) with 
X > 1, — 1 < 2/ < 1 is given by [7] 

ds^ = - f{dt - ujd4>f 



27 / 2 2\( dx^ dy2 



- r -Ft + + - - vn'^r , (1.1) 



where /, u and 7 are functions of x and y only and cr is a constant. They have the 
form 

f - ^„-2qP2Q2 

^ L 

where 



m = ax{l ~ 2;2)(e2«^+ + e^^*- )a+ + yix^ - 1)(1 - a''e''''^^++^-'^)b+ , 

{1 + P2){Qi - Ql) 



1 , , o . x'^ — 1 



7 =-{l + qf\n-^—^ + 2q{l-P2)Qi+q\l-P, 
1 



2 x^ — y 



[X 



2 l)i2Ql ~ 3xQiQ2 + 3Q0Q2 - Q2) 



(1.3) 



2^ 

Here Pi{y) and Qi{x) are Legendre polynomials of the first and second kind 
respectively. Furthermore 

a± - x{l - a2e29(*++5-)) ± (1 + ^2g25(5++5_)) ^ 

b± = ay{e^i^+ + e^^^- ) T ^(6'"*+ - e^^^- ) , 

^± = ^ In + ^(1 - T xy) + J[x(l - y^) T y(-^ - 1)] In , (1.4) 

the quantity a being a constant 
a- M 



a 



ct^Vm^^^. (1.5) 



a 

The Geroch- Hansen [5J [S] moments are given by 

M2fe+i = J2fe-0, A: = 0,1, 2,... (1.6) 

Mo = M, M2^-Ma^ + ^qM^(^l~-^^ , ... (1.7) 

^2x3/2 



Ji=Ma, J3 = _Ma^ + — gM^al 1- — ) (1.8) 
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The vanishing of the odd gravitoelectric (-/l/„) and even gravitomagnetic ( J„) multipole 
moments is a consequence of the reflection symmetry of the solution about the 
hyperplane y = 0, which we will refer to as "symmetry" (or equivalently "equatorial") 
plane hereafter. From the above expressions we see that M is the total mass of the 
body, a represents the specific angular momentum, and q is related to the deviation 
from spherical symmetry. All higher multipole moments can be shown to depend only 
on the parameters M, o, and q. 

We limit our analysis here to the case a > 0, i.e. M > a. In the case a = 
the solution reduces to the extreme Kerr spacetime irrespective of the value of g [2] . 
The case a complex, i.e. a > M, requires a different definition of the quadrupole 
parameter in order to have real Geroch-Hansen moments and can be better discussed 
by using Weyl cylindrical coordinates. We will not explore this case here. 

In this work we analyze some geometric and physical properties of the QM 
solution. The limiting cases contained in the general solution suggest that it can be 
used to describe the exterior asymptotically flat gravitational field of a rotating body 
with arbitrary quadrupole moment. This is confirmed by the analysis of the motion 
of particles on the equatorial plane. It turns out that the whole geometric structure of 
the QM spacetime is drastically changed in comparison with Kerr spacetime, leading 
to a number of previously unexplored physical effects strongly modifying the features 
of particle motion, especially near the gravitational source. In fact, the QM solution 
is characterized by a naked singularity at a; = 1, whose existence critically depends 
on the value of the quadrupole parameter q. In the case q ~ (Kerr solution) x — 1 
represents instead an event horizon. This bifurcating behaviour accounts for the above 
mentioned drastic changes with respect to the Kerr metric. 

Due to the very complicated form of the metric most of the analysis will be 
performed numerically. 



2. Limiting cases 

The QM solution reduces to the Kerr spacetime in the limiting case q ^ and 
to the Erez-Rosen spacetime when a 0. Furthermore, it can be shown that the 
general form of the QM solution (see [Appendix A[ ) is equivalent, up to a coordinate 
transformation, to the exterior vacuum Hartle-Thorne solution once linearized to first 
order in the quadrupole parameter and to second order in the rotation parameter. 



2.1. Kerr solution 

For vanishing quadrupole parameter we recover the Kerr solution, with functions 
c2a;2 + (iV-l ^ (cx + l)(l-2/2) 

JK = 7 , -.sn , ,9 n , tOK — 



{cx + 1)2 + cPy"^ ' c^x"^ + dP-y"^ — 1 

1 / c^x'^ + (fy^ - 

2 c2(x2-y2) 



where 



c=— , d= — , c2 + rf2^1, (2.2) 

so that a — (c— l)/d. Transition of this form of Kerr metric to the more familiar one 

associated with Boyer-Lindquist coordinates is accomplished by the map 

r - M , ^ 

X — , y — cos9 , (2.3) 
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so that X — 1 corresponds to the outer horizon r = 
To first order in q the QM solution becomes 

f = fK- 2q[P2Q2fK + A+S+ + AS-] + 0{q^) 



M 



2P2Q: 



c{x^ - y'^) + x{l - y^) 



1 =1K + q 



-\-ad 



c{x 



2(l-P2)Qi 
y2) + 1 - 



In 



(py'2. _ \ 

^2-1 



1' + ^4 



where 



A± =(P{x±y) 



dPy'^ 



[c{xTy) + i? 



-id. 



- y 
S-) 



0{q' 



(2.4) 



[{cx- 
x±y 



1 



{cd2(a; T 2/)^(l ± + c{x^ - 1)(1 T a;?/) 



a;2 + d2y2_l)2 
(a;Ty)[c2(2:2-l) + rf'(l-y')]} 



(2.5) 



This approximate metric could be used to describe the exterior field of an arbitrarily 
rotating mass source with a small quadrupole moment. The lowest Geroch-Hansen 
multipole moments in this case coincide with those of the exact solution as given 
in Eqs. (|1.7p and (|1.8p . Differences will appear in higher moments where all terms 
containing q^ and higher exponents must be neglected. 

It is interesting to mention that in the hmiting case a — > M, the QM metric leads 
to the extreme Kerr black hole solution, regardless of the value of the quadrupole q [5]. 
This can be easily seen both at the level of the multipole moments (|1.7p and (|1.8p and 
by directly computing the limiting metric. The latter determines in turn the limit of 
applicability of the QM solution, since for values in the range a/M > 1 the multipole 
moments and the metric both become complex. As stated in the Introduction we will 
not consider such a situation here. 



2.2. Erez-Rosen solution 

Similarly, for vanishing rotation parameter we recover the Erez-Rosen solution 
[TUl [TTl [T^ . It is a solution of the static Weyl class of solutions (i.e. oj = 0) with 
functions 

/i^i^ - ^e-2^^=«^ JER = 1, (2.6) 
X + 1 

which reduce to 

x-l _ 1 ^ 

X + 1 2 yx"^ — y 

when q — 0, corresponding to the Schwarzschild solution. 
To first nonvanishing order in a/M we find 



fs-^, ls = ^lni^^] (2.7) 



2{x^-l){x+i) + y^^y - " " y^^y + 

-2(x2-2;2)e2.(^.+^-)] (^) ^ + 0[(„/M)^] , 
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2M |l + ^ [(x + y){y l)e''^'- - {x ~ y){y + l)e''^'-] j - 



+ 0[ia/Mf] , 
7 = 7i.« + ^ {l - 2(^5^ [(1 - y'){e""- + e""-) 

+2{x' - }{^y + 0[(a/M)4] . (2.8) 

This approximate solution can be interpreted as a generalization of the Lense-Thirring 
spacetime [131 [Tl] , which is obtained in the limit g — > by retaining terms up to the 
linear order in a/M. Consequently, the approximate solution (j2.8p could be used 
to investigate the exterior gravitational field of slowly rotating deformed bodies. A 
similar approximate solution, accurate to second order in the rotation parameter and 
to first order in the quadrupole moment, was found long ago by Hartle and Thorne 

m. 



2.3. Hartle-Thorne solution 



The Hartle-Thorne metric describing the exterior field of a slowly rotating slightly 
deformed object is given by 

-1 



1 



2M \ 
R ) 



l + 2fciP2(cose) + 2 1 



2M\ 
R ) 



i?4 



f2cos2e - 1' 



1-^^ l-2A:2P2(cose) 
i?2(de2 + sin2 QA(t?)[l - 2fc3P2(cose) 



R ) i?4 

J . 2 



dR' 



4— sin^ Qdtd4> , 
R 



(2.9) 



where 



fci = 



MR^ 



M 
R 



1 + ^ 



5 Q - J^/M , 



k2 — ki — 



k3 = ki + ^-- 



6J^ 

5 Q - jyM 



M^R 



2M 
~R 



-1/2 



1 - ^ Q\ 



M 



.(2.10) 



Here are the associated Legendre functions of the second kind and the constants 
M., J and Q are the total mass, angular momentum and mass quadrupole moment 
of the rotating star respectively. 

It is interesting to find out the connection between the Hartle-Thorne solution 
and the QM solution in the appropriate limit. To this end it is necessary to start with 
the general form of the QM solution as given in [Appendix A| containing an additional 
parameter, the Zipoy-Voorhees [THl E] constant parameter 6. For our purposes it 
is convenient to set such a parameter as (5 = 1 -I- sq, where s is a real number. For 
5 = 0, i.e. (5 = 1, we recover the solution (|l.ip - (|1.5p . To first order in the quadrupole 
parameter q and to second order in the rotation parameter ajM the metric functions 
turn out to be 



/ 



1 



\-q \ 2P2Q2 - sin- 



1 



2y^ 



{x + lf 



f-V 



6 



uj ~ 2M 



1 \m) 



II -y 



7 ~7(l + 2s<7)- -^2 



1 \m) ' 



(2.11) 



where 7 is defined in Eq. (|1.3p and terms of the order of q{a/M) have also been 
neglected. 

Introduce first Boyer-Lindquist coordinates {t, r, 6, (p) through the transformation 

r-M 

X — , y = cos , (p — (p . 



Then set s = — 1 and M = M'{1 + q). The further transformation 

t^t , r = r{R,e) , e = 9{R,e) , 

with 



(2.12) 
(2.13) 



r = R + M'q + -M' q sin^ 9 



R 1 i?2 

M' ^ ^ 2M^ 



1 - 



2M' 



2M' 



a 

2^ 



= 8 — sin 8 cos 8 < -o 
' 2 



1-^1- cos2 8(1- 



R 



R 
2M' 



In 1 - 



2M' 



— - 1 1 In ( 1 - 
M' 



R 
2M' 



1 



i? 



a 

2^ 



1 



2M' 



(2.14) 



finally gives the mapping between the general form of QM solution and the Hartle- 
Thorne metric (j2.9p with parameters 

M - M' = M(l - g) , = -A/a , Q=^ + -M^q . (2.15) 

Al 5 

Note that the previous transformation is obtained simply by combining the 
corresponding transformation from Kerr to Hartle-Thorne solution as given by Hartle 
and Thorne themselves in [15 and that from Erez-Rosen to Hartle-Thorne solution 
as found by Mashhoon and Theiss in [18]. 



3. Geometric properties of the solution 

The solution admits two Killing vectors associated with translation in time and 
rotation about the symmetry axis. The timelike Killing vector dt changes its causality 
property when / = 0, defining a hypersurface which in the Kerr limiting case is called 
ergosurface as the boundary of the ergosphere (or equivalently ergoregion). Using 
this terminology also in this case we compare in Figs. [1] and [5] the ergoregions of 
Kerr solution and QM solution (for different values of the quadrupole parameter) 
respectively. The situation is completely different. In fact, while in the Kerr case such 
a hypersurface is the boundary of a simply-connected domain, in the case of the QM 
solution this property is no more true, as soon as the magnitude of the quadrupole 
parameter exceeds a certain critical value (for example, for the choice a/M — 0.5 the 
range of Kerr-like behaviour corresponds to — 1 < q < 1.21). 

The spacelike Killing vector also changes its causality condition in some 
regions, leading to the existence of closed timelike curves in the QM solution. These 
regions are depicted in Fig. [3] for different values of the quadrupole parameter. 

In order to investigate the structure of singularities in the QM metric we have to 
consider the curvature invariants. Since the solution is a vacuum one there exist just 
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two independent quadratic scalar invariants, the Kretschmann invariant A'l and the 
Chern-Pontryagin invariant K2 defined by 

= R'^P^^R^p-.s , K2 = *R''^''^Ro,0^5 , (3.1) 

where the star denotes the dual. The behaviour of curvature invariants as functions 
of X for selected values of y, a/M and q is shown in Fig. 01 It turns out that Ki 
diverges when approaching x = 1 along the direction y = 0, but it is finite there 
moving along a different path. Furthermore, the invariant K2 vanishes identically for 
y = 0. Therefore, the metric has a directional singularity at a; = 1 (see e.g. [H]). 

It is now interesting to investigate what kind of hypersurface is a; = 1. Consider 
the normal to a a; = const hypersurface. The behaviour of its norm g^^ when 
approaching x = 1 discriminates between its character being either null or timelike 
depending on the value of the quadrupole parameter. The limit of as a; — > 1 also 
depends on y. For instance, approaching a; = 1 along any direction on the equatorial 
plane y = gives 

- (x- l)i+«/2-9V4 ^ (3 2) 

implying that the singular hypersurface a; = 1 is null when \q— 1| < \/5 and timelike 
otherwise. On the other hand, moving along the axis y — 1 gives 

r (a: - 1)1+9 , g > 

g-^; -l<g<0 (3.3) 

I {x-ir'i , g<-l , 

implying that the singular hypersurface a: = 1 is always null. The corresponding 
careful analysis for the Erez-Rosen solution (a = 0) is contained in Appendix A of 

120]. 

A similar discussion can be done also for the metric determinant. Approaching 
a; = 1 along any direction on the equatorial plane y = gives 

,/^^(a:-l)-9/2+9^/S (3.4) 

implying that the volume element vanishes approaching the singular hypersurface 
a: = 1 when q < 0, q > 2 and diverges otherwise. On the other hand, moving along 
the axis y — 1 gives 

r ix-l)-'^ , q>0 
^ < const , -1 < q < (3.5) 

[ (a; - l)«+i , q<-l . 

We note that if the presence of the quadrupole moment totally changes the situation 
with respect to the Kerr spacetime, the smooth horizon of the Kerr solution becoming 
a singular hypersurface in the QM solution, the properties of the naked singularity 
are also different with respect to the limiting case of the Erez-Rosen spacetime. 

Finally, the spectral index S (see [Appendix B| for the definition of S in terms 
of Weyl scalars) shows that the solution is algebraically general. The real part of 
S is plotted in Fig. [5] as a function of x for ?/ = and selected values of a/M and 
q. The imaginary part is identically zero in this case. A numerical analysis of the 
spectral index for different values of the quadrupole moment shows that S" ^ 1 as 
X ^ Qo, i.e. it is algebraically special. We conclude that the asymptotic behavior 
of the spacetime is dominated by the Kerr spacetime which is algebraically special 
of type D. This is in agreement with the expectation which follows from the analysis 
of relativistic multipole moments, according to which any stationary axisymmetric 
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.2 1.4 1.6 18 

X 



Figure 1. The shape of the ergoregion is shown in the case of vanishing 
quadrupole parameter q = (Kerr spacetime) for the choice a/M = 0.5 in prolate 
spheroidal coordinates. 



asymptotically flat vacuum solution of Einstein's equations must approach the Kerr 
metric asymptotically 

Note that drawing a Penrose diagram would help to easier understand the global 
aspects of the geometry of the QM solution. However, due to the rather involved form 
of the metric functions it is a very hard task to construct it. Analytic calculations can 
be performed only in the case of small values of the quadrupole parameter q. But in 
this simplest case ((/ <C 1) it is possible to show that the resulting conformal diagram is 
closely similar to the corresponding one for a Kerr naked singularity (see e.g. Fig. 6.4 
ofRef. [21] )• 



4. Geodesies 



The geodesic motion of test particles is governed by the following equations: 



■ _ E ujf 



(L - LoE) , 



f 



1 y2 

2X2 



+ 



1 



fy 
f 



27y 



x2 + r2 

X2 + y2 y2 



/ ^^"^ X2 + y2 



1 



-27 



2 fa'^X^Y^iX^ +Y^) 
+E'^a'^X^Y^]fy + 2{L - ujE)f[y{L - luE) - EY^Uy]] , 



-27^2 



X' 



y2" cr2(X2 + r2) 



E^ - fl\f - 



P 



2^2y: 



;{L-UjEf 



xy 

{Y'[f\L~ojEf 
(4.1) 



where Killing symmetries and the normalization condition have been used. Here E 
and L are the energy and angular momentum of the test particle respectively, ^ is 
the particle mass and dot denotes differentiation with respect to the affine parameter; 




(c) (d) 

Figure 2. The shape of the ergoregion is shown for ajM = 0.5 and different values 
of the quadrupole parameter: q = [—10,-1,1,10], from (a) to (d) respectively. 
For q = —1 and q = I the shape is similar to the Kerr case, i.e. the regions are 
simply-connected. 
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(a) (b) 

Figure 3. The regions wtiere ttie metric component g^^ changes its sign are 
shown for a/M = 0.5 and different values of the quadrupole parameter: (a) 
q = 10 and (b) q = —10. The existence of closed timclikc curves is allowed there. 



furthermore, the notation 



X = yjx"^ - 1 , Y = v/1 (4.2) 

has been introduced. 

Let us consider the motion on the symmetry plane y — Q. If y = and y = 
initially, Eq. (j4.ip ?{ ensures that the motion will be confined on the symmetry plane, 
since /y, ujy and 7^ all vanish at y = 0, so that y = too. Eqs. (14. ip thus reduce to 



.2 6-27^2 



fT2(l+X2) 



E'-f^'f-J^iL-cuEf 



(4.3) 



where metric functions are meant to be evaluated at y = 0. The motion turns out to 
be governed by the effective potential V defined by the equation 

- - ^ {L-u;Vf^O. (4.4) 



In fact, for E — V the rhs of Eq. (|4.3P q vanishes. 

The behaviour of V as a function of x is shown in Fig. [51 Repulsive effects occur 
for decreasing values of x approaching x = 1. It is worth to mention that the repulsive 
effect of the naked singularity has also been investigated by Herrera [H] in the case of 
a quasispherical spacetime but without taking into account the rotation of the source. 

The case of a geodesic particle at rest will be analyzed below. Circular geodesies 
will be discussed in detail in the next section, where accelerated orbits are also studied. 
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(c) 

Figure 4. The behaviour of the Kretschmaiiii invariant Ki as a function of x 
is shown for tlie choice of parameters a/M = 0.5 and q = 10 for different values 
of y: (a) y = and (b) y = 0.5. Figure (c) shows instead tlie beliaviour of the 
Chern-Pontryagin invariant K2 for the same clioicc of parameters as in Fig. (b). 
Figs, (b) and (c) show how both the invariants Ki and K2 change their signs as 
X approaches unity. This is a manifestation of the appeaxence of repulsive gravity 
regions, typical of naked singularity solutions. When a; — > 00 the invariants Ki 
and K2 both vanish, as expected. 
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Re(S) 3 




2 3 4 5 6 

X 



Figure 5. The behaviour of the real part of the speciahty index S as a function 
of X is shown for the choice of parameters a/A/ = 0.5, g = 10 and y = 0. In this 
case (j/ = 0) we have also Im(S') = 0. 



^.1. Particle at rest 

For the QM solution which is characterized by the presence of a naked singularity it is 
even possible to satisfy the conditions for a geodesic particle to be at rest for certain 
values of the quadrupole moment. This effect can be explained when the attractive 
behaviour of gravity is balanced by a repulsive force exerted by the naked singularity. 
Usually, repulsive effects are interpreted as a consequence of the presence of an effective 
mass which varies with distance and can thus become negative. The consideration of 
the corresponding post-Newtonian limit shows that an effective mass can be indeed 
introduced, depending on the distance from the source and the value of the Geroch- 
Hansen quadrupole moment [23| . We see that in the case of the corresponding exact 
solution under consideration a similar situation takes place. 
A particle at rest is characterized by the four velocity 

U - -^dt . (4.5) 

The corresponding four acceleration a{U) — VjjU is given by 

= 2^2(^"2 + y2) [-^'/^^^ + Y'fydy] • (4.6) 

On the symmetry plane y — we have fy — 0, so that the geodesic condition 
a{U) = implies — 0. The pairs {x,q) satisfying this condition are shown in 
Fig. [7] for a fixed value of a/M. As an example, for a/M = 0.5 and g = 10 we 
get X « 1.588 as the equilibrium position. In this case the corresponding energy 
and angular momentum per unit mass of the particle are given by E/ fi « 0.553 and 
L/{iiM) « 0.257 respectively. No equilibrium positions exist for q<0. 
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Figure 6. The behaviour of the effective potential V as a function of x is shown 
for the choice of parameters a/M = 0.5 and L/{fj,M) = 10 for different values of 
the quadrupole parameter: (a) q = 0, (h) q = 10 and (c) q = —10. 
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(a) 



(b) 



Figure 7. The pairs (x, q) allowing a test particle to be at rest in the QM 
spacotime are shown for a/M = 0.5. Fig. (b) is a detail of Fig. (a) close to a; = 1. 
No equilibrium positions exist for q < 0. 



5. Circular orbits on the symmetry plane 

Let us introduce the ZAMO family of fiducial observers, with four velocity 

n = N-Hdt-N'^d^); (5.1) 
here N = (-5**)-^/^ and iV^ = gt,p I g^fj, are the lapse and shift functions respectively. 



A suitable orthonormal frame adapted to ZAMOs is given by 



1 



— ^1 



--dx 



1 



1 



(5.2) 



/ 9xx y 9yv \ 

with dual 

J = NAt , Lu^ = Vft^dx , ujy = V5^dy , J = y^idq^ + N^dt) . (5.3) 

The 4- velocity U of uniformly rotating circular orbits can be parametrized either 
by the (constant) angular velocity with respect to infinity ^ or, equivalently, by the 
(constant) linear velocity i' with respect to ZAMOs 

U^m+Cd4=^[ei + i^e^], 7 = (l-z/2)-i/2 ^ (5 4) 

where F is a normalization factor which assures that UqU°' = — 1 given by: 
and 

N 



(5.5) 



C = -N'f' 



/9<t><P 



(5.6) 



We limit our analysis to the motion on the symmetry plane (jj = 0) of the solution 
(|l.ip - p.5p . Note both y = and x — xq are constants along any given circular orbit. 
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and that the azimuthal coordinate along the orbit depends on the coordinate time t 
or proper time t along that orbit according to 

- <^o = a = ^uru , ^u^TC, (5.7) 

defining the corresponding coordinate and proper time orbital angular velocities C and 
flu- These determine the rotation of the spherical frame with respect to a nonrotating 
frame at infinity. 

The spacetime Frenet-Serret frame along a single timelike test particle worldline 
with 4-velocity U = Eq and parametrized by the proper time tu is described by the 
following system of evolution equations [H] 

DEo DEi 

— — = kEi , — — = kEq + T1E2 , 

DE2 DE3 , , 

— ^ = -ri£;i + T^E^ , — ^ = -r^E^ . (5.8) 
drij dTu 

The absolute value of the curvature k is the magnitude of the acceleration a(C/) = 
DU/dru = K-Ei, while the first and second torsions ri and T2 are the components of 
the Frenet-Serret angular velocity vector 

^(PS) - riS3 + r^E^ , ||o^(FS)ll = [rl + r^Y'^ , (5.9) 

with which the spatial Frenet-Serret frame {Ea\ rotates with respect to a Fermi- 
Walker transported frame along U . It is well known that any circular orbit on the 
symmetry plane of a reflection symmetric spacetime has zero second torsion T2 , while 
the geodesic curvature n and the first torsion ti are simply related by 

1 d/i , , 

On the symmetry plane there exists a large variety of special circular orbits 
[25l [26l [27l [28j : particular interest is devoted to the co-rotating (+) and counter- 
rotating (— ) timelike circular geodesies whose linear velocities is 

''te™)* ^± V^^a{/,[/2a;2 + a2(a;2-l)]+2/(/W-a2a;)} ' 
where 

C = - 2(72(a;2 - l)Lof., - /{tu^lfuj^ + C7^{x'' - 1)] - 2a^xuj} , 

D = fLol - a^MM^^ - 1) - 2x/] . (5.12) 

All quantities in the previous expressions are meant to be evaluated at y = 0. The 
corresponding timelike conditions \h'±\ < 1 together with the reality condition D > 
identify the allowed regions for the "radial" coordinate where co/countcr-rotating 
geodesies exist. 

Other special orbits correspond to the "geodesic meeting point observers" defined 
in [27], with 

^(gmp) = 2 ■ ^ ' 

A Frenet-Serret (FS) intrinsic frame along U [24] is given by 
Eq = U ^ "fin + fe^] , Ei = e±, E2 ^ Cy , E2 = E, ^ ^[vn + e^] . (5.14) 
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It is also convenient to introduce the Lie relative curvature of each orbit [57] 
^(lic) = - 9i In y^g^ 

2aa;V7 I a^x^ - 1) - poj^ 

It then results 

K = fc(iio)7^('^ - '^+){'^ - >^-) , 

n = fc(lio)J^(gmp)7^('^ - i^(crit) + )('^ - i^(crit)-) , (5.16) 

where 

_ 1-1^- T 1+1^+ . 

^(crit)± — — , ^^(crit) + i^(crit)- — ^ ) (^■^' 

7- T 7+ 

identify the so called "extremely accelerated observers" [371 US]: i'(ext) = ^^(crit)-) 
which satisfies the timelike condition in the regions where timelike geodesies exist, 
while J^(crit)+ is always spacelike there. 

The geodesic velocities (|5.1ip are plotted in Fig. [8] both as functions of the 
quadrupole parameter q for fixed "radial" distance (see Fig. (a)) and as functions 
of X for different values of q (see Fig. (b)). In the first case (Fig. (a)) we have shown 
how the quadrupole moment affects the causality condition: there exist a finite range 
of values of q wherein timelike circular geodesies are allowed: qi < q < qg, for co- 
rotating and q2 < q < q3 for counter-rotating circular geodesies. The critical values 
qi , 52 and q^ of the quadrupole parameter can be (numerically) determined from Eq. 
(|5.1ip . The difference from the Kerr case is clear instead from Fig. (b): the behaviour 
of the velocities differs significantly at small distances from the source, whereas it is 
quite similar for large distances. 

A similar discussion concerning the linear velocity '^(gmp) of "geodesic meeting 
point observers" as well as J^(ext) of "extremely accelerated observers" is done in Figs.|9] 
and do] respectively. 

The behaviours of both the acceleration k and the first torsion ti as functions 
of ly are shown in Figs. [Tl] and [T2] for different values of the quadrupole parameter 
and fixed x as well as for different values of the "radial" distance and fixed q. A 
number of interesting effects do occur. For instance, from Fig. [Tl] (a) and (b) one 
recognizes certain counter-intuitive behaviours of the acceleration, in comparison with 
our Newtonian experience. These effects have their roots also in the Kerr solution and 
are well known and studied since the 90s ^27j. For instance, for negative values of q 
{q = —500, —250) we see that by increasing the speed (for positive values over that 
corresponding to the local minimum) the acceleration also increases; hence, in order 
to maintain the orbit, the particle (a rocket, say) should accelerate outwards. This 
is counter-intuitive in the sense that for a circular orbit at a fixed radius increasing 
the speed corresponds to an increase of the centrifugal acceleration and therefore to 
maintain the orbit "classically" one would expect to supply an acceleration inward. All 
such effects have been analyzed in the past decades in the Kerr spacetime in function 
of the radius of the orbit, i.e. the distance from the black hole. The novelty here 
is represented by Fig. [TT] (a) where the various curves do not correspond to different 
orbital radii (i.e. different values of the coordinate x, as it is for the cases (b) and 
(c)), but to different values of the quadrupole parameter g at a fixed radial distance. 
Therefore, the conclusion is that a spacecraft orbiting around an extended body - 
according to general relativity and the QM solution- should expect counter-intuitive 
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(a) (b) 

Figure 8. The geodesic linear velocities '^^^^-^ a-re plotted in Fig. (a) as functions 
of the quadrupole parameter q for fixed distance parameter x = 4 from the source 
and a/M = 0.5. Co-rotating and counter-rotating circular geodesies exist for 
qi < q < qs and 92 < 9 < <?3 respectively, with qi ^ —105.59, 92 ~ —36.29 and 
93 ~ 87.68 for this choice of parameters. The behaviour of f (^g^) as functions of 
X is shown in Fig. (b) for different values of q = [-80, -30, -10, 0, 2, 4, 10, 50]. 
The thick curves correspond to the Kerr case (</ = 0). Curves corresponding to 
great positive value of the quadrupole parameter in the allowed range exhibit both 
a local maximum (I'^^g^j) and a local mimimum i^^^^^^)- For decreasing values 
of q the local minimum first disappears; for q further decreasing also the local 
maximum disappears (the curves are thus ordered from left to right for increasing 
values of q). Curves corresponding to negative values of q never present extrema 
as in the case of Kerr spacetime (g = 0), the lightlike condition being reached at 
greater values of the "radial" distance for decreasing values of q. 



engine acceleration to remain on a given orbit. Outward or inward extra acceleration 
for an increase of the speed critically depend on the quadrupole moment (i.e. the 
physical structure) of the source, a fact that should be taken into account. 

Fig. [13 shows instead that at a fixed radius one can always find a value of the 
quadrupole parameter and a value of the speed at which the first torsion vanishes. 
Being the second torsion identically vanishing, in these conditions the Frenet-Serret 
frame becomes also a Fermi- Walker frame: in fact the Frenet-Serret angular velocity 
represent the rate of rotation of the Frenet-Serret frame with respect to a Fermi- 
Walker one. In (see Sec. 6, Eq. (6.22)) the precession of a test gyroscope 
(Fermi- Walker dragged along a circular orbit) has been related to the first torsion of 
an observer-adapted Frenet-Serret frame. The same discussion repeated here together 
with a simple inspection of Fig. [12] allows us to include the effects of the presence of 
the quadrupole parameter q. 
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(a) (b) 

Figure 9. The linear velocity i^(gnip) corresponding to the "geodesic meeting 
point observers" is plotted in Fig. (a) as a function of the quadrupole parameter 
q for fixed distance parameter x = A from the source and a/M = 0.5. The 
behaviour of t'(ginp) ^ ^ function of x is shown in Fig. (b) for different values 
of g = [—50, —10, 0, 2, 5, 10, 30, 80]. It exists only in those ranges of x where both 
'^(gco) '^'^'^ '^(gco) ^^^^^ (solid black lines). These ranges are listed below; q = —50, 
X > 4.23; q = -10, x > 3.36; q = 0, x > 2.92; g = 2, 1.0012 < x < 1.13 and 
X > 2.79; 5 = 5, 1.21 < x < 1.64 and x > 2A8; q = 10, x > 1.54; q^= 30^ x > 2.44; 
g = 80, a; > 3.83. ~ ~ ~ ~ ~ 

6. Concluding remarks 

In this work we investigated some properties of the QM spacetime which is a 
generahzation of Kerr spacetime, including an arbitrary mass quadrupole moment. 
Our results show that a deviation from spherical symmetry, corresponding to a 
non-zero gravitoelectric quadrupole moment, completely changes the structure of 
spacetime. A similar behaviour has been found in the case of the Erez-Rosen 
spacetime [30] . A naked singularity appears that affects the ergosphere and introduces 
regions where closed timelike curves are allowed. Whereas in the Kerr spacetime the 
ergosphere corresponds to the boundary of a simply-connected region of spacetime, in 
the present case the ergosphere is distorted by the presence of the quadrupole and can 
even become transformed into multiply-connected regions. All these changes occur 
near the naked singularity. 

The presence of a naked singularity leads to interesting consequences in the motion 
of test particles. For instance, repulsive effects can take place in a region very close 
to the naked singularity. In that region stable circular orbits can exist. The limiting 
case of static particle is also allowed, due to the balance of the gravitational attraction 
and the repulsive force exterted by the naked singularity. 

We have studied the family of circular orbits on the symmetry plane of the 



19 




(a) (b) 

Figure 10. The linear velocity V((,xt) corresponding to the "extremely accelerated 
observers" is plotted in Fig. (a) as a function of the quadrupole parameter 
q for fixed distance parameter x = A from the source and a/M = 0.5. The 
behaviour of fjext) ^ ^ function of x is shown in Fig. (b) for different values of 
q = [—50,-10,0,5,10,30,80]. The dashed curve corresponds to the case q = 0. 
It exists only in those ranges of x where both i'^ and 1^7^^^ exist (see Fig. [Q}. 



QM solution, analyzing all their relevant intrinsic properties, namely Frenet-Serret 
curvature and torsions. We have also selected certain special circular orbits, like the 
"geodesic meeting points" orbits (i.e. orbits which contain the meeting points of two 
oppositely rotating circular geodesies) and the "extremely accelerated" orbits (i.e. 
orbits with respect to which the relative velocities of two oppositely rotating circular 
geodesies are opposite) , whose kinematical characterization was given in the 90s [27] 
with special attention to Kerr spacetime. Here we have enriched their properties 
specifying the dependence on the quadrupole parameter. 

The question about the stability of the QM solution is important for astrophysical 
purposes. In this context, we have obtained some preliminary results by using 
the variational formulation of the perturbation problem as developed explicitly by 
Chandrasekhar |31j for stationary axisymmetric solutions. A numerical analysis 
performed for fixed values of the parameters entering the QM metric shows that 
it is unstable against perturbations that preserve axial symmetry. One can indeed 
expect that, once an instability sets in, the final state of gravitational collapse will be 
described by the Kerr spacetime, the multipole moments of the initial configuration 
decaying during the black hole formation. Nevertheless, a more detailed analysis is 
needed in order to completely estabhsh the stability properties of this solution. 

Finally, we mention the fact that it is possible to generalize the metric investigated 
in this work to include the case of a non spherically symmetric mass distribution 
endowed with an electromagnetic field [32] . The resulting exact solution of Einstein- 
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(a) (b) 
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Figure 11. The acceleration re for circular orbits at 1/ = is plotted in Fig. (a) 
as a function of v for ajM = 0.5, x = A and different values of the quadrupole 
parameter: q = [—.500, —250, —100, 0, 80, 250, 500]. The curves are ordered from 
top to bottom for increasing values of q. The values of u associated with k = 
correspond to geodesies, i.e. ^(^eo)' "^^^ behaviour of k as a function of v 
for different x is shown in Figs, (b) and (c) for fixed values of the quadrupole 
parameter: (b) g = 1, x = [1.1, 1.25, 1.5,2.5,4] and (c) q = 100, x = [2.5,3,4, 10]. 

Maxwell equations turns out to be asymptotically flat, contains the Kerr-Newman 

black hole spacetime as a special case, and is characterized by two infinite sets of 
gravitational and electromagnetic multipole moments. For a particular choice of the 
parameters, the solution is characterized by the presence of a naked singularity. It 
would be interesting to explore repulsive effects generated by the electromagnetic field 
of the naked singularity also in this case. This task will be treated in a future work. 
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Figure 12. The first torsion t\ for circular orbits at j/ = is plotted as a function 
of V for a/M = 0.5, x = A and different values of the quadrupole parameter: (a) 
q = [0,100,250,500,750], (b) q = [-1000,-750,-500,-250,-150] and (c) q = 
[—150, —80, 0]. The curves are ordered from left to right for increasing/decreasing 
values of q in Fig. (a)/Fig. (b) respectively. Figure (c) shows the changes of 
behavior occurring in the different regions <? < gi, <?i < q < q2 and q > q2 (see 
Fig. H. 
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Appendix A. General form of QM solution with arbitrary Zipoy-Voorhees 
parameter 

The general form of QM solution with arbitrary Zipoy-Voorhees parameter S is given 
by Eq. with functions 

f _ ^p-2qSP2Q2 

L 



m 

u> = —2a — 2a— e 
R 

4 V a 



2qSP2Q2 

R 



If. M I 

(a;2 - 1)^ 

where 7 is the same as in Eq. (|1.3p . while 

R = a^a^ + , 



L = al + bl , 

Tl={x + l)^~'^[x{l - y^){X + ^l)a+ + y{x^ - 1)(1 - \^l)b^ 
The functions a± and b± are now given by 

a± = {x± 1)''"^ - A^i) ± (1 + A^)] , 
6± ={x±lf-\x{\ + ^Ji)T[\-^J)] , 

with 



\25-2 2qSS^ 



(A.l) 



(A.2) 



(A.3) 



A ^ a{x' - ly-'ix + y) 

/i = a{x^ - iy-'{x - y)^S-2^2qSS^ (A 4) 

The functions 6± and the constants a and a arc instead the same as in Eqs. (|1.4p and 
p.Sp respectively. 

This solution reduces to the solution (|l.ip ~ (|1.5p for 5 = 1. 



Appendix B. Newman- Penrose quantities 

Let us adopt here the metric signature (+,—,—,—) in order to use the Newman- 
Penrose formalism in its original form and then easily get the necessary physical 
quantities [3T] . The Weyl-Lewis-Papapetrou metric is thus given by 

ds^ = f{dt - Ujd(j> f 



dx'^ 



Introduce the following tetrad 
I 



dy^ 



1 l-y^ 



+ {x^-l){l-y^)dcl,' 



n = \ - 



where 





dt 








dt 






1 a 




\/ x^ 


- 1 



aXY\ 

+ — ) 

aXY\ 



Y2 



dx .dy 



(B.l) 



(B.2) 
(B.3) 
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The nonvanishing spin coefficients are 
K = ~ A 

T = -TT* = A{xY - iyX) , 
V = - A 



XY f 

— (X/, + iVfy) + ^{Xco^ + iViVy) +xY- iyX 
J 



XY f 

— (X/^ - iYfy) + ^{Xu^ - iYwy) -xY- iyX 
J ^ 



a = AXYa'^ 
+ 



2/ 

xX + iyY 



UXU - iY!y) + - lY^y -£{y- i^) 



x2 + y2 

P = AXYa^ ^{Xf, + zYfy) - - lYjy - ^ + » y) 
xX — iyY 



X^ + Y^ 



where 



A^iV2-^'~' 



I 



4 (jXY \ X^ + Y"^ ' 
The nonvanishing Weyl scalars are 



(B.4) 



(B.5) 



ipo = AXY 



t{XU + iYfy) + 2f{XK^ + iYKy) 



aXY 



{k + t){Xu)x + iYujy) 



+ 2k/3 + t{k - t) , 
1 A^X^Y^ 



2 = 7: 



x[ fx + 2 fix - 3i 



aX^ 



y{fy + Vly + 3? 



6 / 
f 

+ 3i-{fxUJy - fyUx) - X^{f^^ - 2/7^^) - Y'^{fyy - 2f^yy) 

a 

C72 I y2 ^ X2 



P 



(7^2 



-7r(X/, - iYfy) - 2f{Xu, - iYvy) - -^(y + 7r)(Xa;, - iYu^y) 



= AXY 

+ 2i/a + TT{u - tt) . (B.6) 

Finally the two scalar invariants of the Weyl tensor whose ratio defines the 
speciality index have the following expressions in terms of the Newman-Penrose 
curvature quantities 



I = VoV'4 + 3V'2 , J = ■>Poi)2iPi -n, S = 27 

S has the value 1 for algebraically special spacetimes. 



(B.7) 
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